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Abstract. - Unlike in ordinary metals, in graphene, phonon structure can be seen in the quasi- 
particle electronic density of states, because the latter varies on the scale of the phonon energy. 
In a magnetic field, quantization into Landau levels creates even more significant variations. We 
calculate the density of states incorporating electron-phonon coupling in this case and find that 
the coupling has pronounced new effects: shifting and broadening of Landau levels, creation of 
new peaks, and splitting of any Landau levels falling near one of the new peaks. Comparing our 
calculations with a recent experiment, we find evidence for a phonon with energy similar to but 
somewhat greater than the optical E2 g mode and a coupling corresponding to a mass enhancement 
parameter A ~ 0.07. 



Nature has presented us with a fascinating solid state 
system in the form of graphene. This material is a two- 
dimensional crystal of carbon atoms arranged on a honey- 
comb lattice, an apparently simple structure that gives rise 
to astonishing electronic and mechanical properties [l}[3]- 
The charge carriers display a linear dispersion at low ener- 
gies, with the low-energy limit of the tight-binding Hamil- 
tonian mapping onto a relativistic Dirac Hamiltonian for 
massless fermions, specifically the Weyl Hamiltonian for 
massless neutrinos |ij[5j . In this Hamiltonian, the valence 
and conduction bands at low energies intersect linearly at 
a point called the Dirac point (DP). Angle-resolved phote- 
mission spectroscopy (ARPES) experiments have observed 
this band structure [6] , and correspondingly, scanning tun- 
neling spectroscopy (STS) experiments have measured an 
electronic density of states (DOS) that grows linearly out 
of the DP [Tj [s] . However, there is clear evidence that 
this behavior is altered by many-body renormalizations: 
in both the ARPES and STS data, the slopes of the lin- 
ear curves are modified, and kinks or bumps exist at fre- 
quencies corresponding to phonon modes [6{ ff\ [9}fl2] or 
to electron-hole excitations or electron-plasmon interac- 
Evidcnce of many-body renormalizations is 
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tions 

also seen in other experiments, such as optical conduc- 



bare predictio n [14||15] , attributed to electron-electron in 
teractions 16p7 and electron-phonon 18p9 
that create Holstcin sidebands; see ref. 20 



interactions 
for a review. 



In this letter, we focus on one particular many-body 
effect: phonon signatures in the electronic DOS. In ordi- 
nary metals, these signatures are not seen at all, because 
the DOS is essentially energy-independent. Conversely, 
in graphene, the DOS varies significantly on the scale of 
the phonon energy, meaning phonon structure can be seen 
in it 19 21 . This has allowed scanning tunneling spec- 



troscopy experiments to measure bumps and changes in 
slope of the DOS (or its first derivative) at energies typi- 
cal of phonons in graphene, specifically the acoustic out- 
of-plane mode at 67meV 
around 165meV 
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tivity measurements 13 that find absorptions beyond the 



and the Ai g optical mode 
(The E2 g optical mode around 
200meV, which has been seen in ARPES data [6j, has not 
been observed in STS, though a phonon structure of un- 
known origin has been observed at the somewhat higher 
value of 240meV 12 ; this will be discussed below in con- 
nection with our own results.) However, evidence still 
remains somewhat limited. In a magnetic field 22-24 
there is much more opportunity to determine the effects 
of electron-phonon coupling, because the DOS condenses 
into Landau levels (LLs), making the energy-dependence 
of the DOS much more extreme. Moreover, the levels can 
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Fig. 1: (Color online) (a) Bare (dashed red curve) and renormalized (solid black) density of states and self-energy as a function 
of energy, with A cff = 0.04, loe = 200meV, vf = 10 6 m/s, T = 5meV, fi = 0, and B = 15T. For this case of zero chemical 
potential, the negative-energy region is given by N(—lj) = N(w), Ei(— u>) = — Ei(u>), and E2(— ui) = E2(w). The DOS is in 
units of NoeV, and E is in units of meV. (b) DOS for increasing values of B, with parameters otherwise the same as at left. 
Black arrows point to the fourth Landau level, which is split in two when crossing the phonon energy. Curves are offset for 
clarity, (c) A schematic representation of the position of bare (left) and renormalized (right) energy levels on the Dirac cone. 
Landau levels are shown as numbered disks outlined in dark blue. Phonon peaks are indicated by the thick dashed red curves. 
Any Landau level sufficiently close to a phonon peak will be split in two: in this case, n = ±2 — > ±2 . 



be moved by adjusting the strength of the magnetic field. 
In an ordinary two-dimensional electron gas, the LLs are 
evenly spaced by the cyclotron frequency uj c — eB/mc, 
but in graphene their energies are proportional to y/\n\B, 
where n is an integer. In addition, in graphene the n = 
LL is field-independent and fixed to the DP, where there 
are equal parts hole and electron states [4j |T] . These 
LLs have been clearly observed in experiments. In par- 
ticular, tunneling measurements have confirmed the field- 
dependence [7,8:25 27] , and optical conductivity measure- 



ments [28] have observed the predicted 23[ 24 
between LLs. 



transitions 



Because the magnetic field can moderate the DOS con- 
siderably through the LLs, the amplitude of phonon struc- 
tures should be highly dependent on B. Moreover, the 
quantization of the LLs in graphene in typical magnetic 
fields used in experiment is on the scale of the phonon 
energy, making it possible to sweep the LLs through the 
phonon energy by changing B. All of these features have 
great potential to identify phonon structures, the energy 
of the phonon, and the electron-phonon mass renormal- 
ization parameter A. In this letter, we characterize the 
signatures of phonons in the DOS of graphene in a mag- 
netic field: renormalization of LL energies, splitting of 
LLs, phonon structures, and other signatures. We fur- 
ther investigate these ideas by providing a comparison of 
the theory with recent data of Miller et al. 18] . 

In the absence of a magnetic field, the electronic density 
of states in graphene varies linearly as N(co) = Nq\uj + /zo|, 
where No = -¥§-a and /io is the non-interacting chemi- 
cal potential. Once a magnetic field B is applied, and 
accounting for many-body renormalizations due to an 



electron-phonon self-energy £ = + i^2i the density of 
states becomes a sum of Lorentzians, one for each Landau 
level 
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[u - Ei(w) + n - M„] 2 + [r - S 2 (w)] 2 ' 



(1) 



where M n = sgn(n) yj2\n\eBv 2 F h/ c is the energy of the nth 
Landau level, /i — /iq + Xi(0) is the interacting chemical 
potential [19' , is a high-energy cutoff (taken to be 7eV 
throughout this letter) , and T is a residual scattering rate 
(taken to be constant for simplicity). At zero temperature 
the self-energy is given by |19][2l] 



E(«) 



TF 



dva 2 F(v) \ duj 



,7V(o/) 
0(-u/ 

v — U)' 



(2) 



where a 2 F(y) is the electron-phonon spectral density. For 
the moment we assume an Einstein phonon mode at loe, 
described by a 2 F(v) — A8{y — lje), where A is the 
electron-phonon coupling strength. Such a phonon spec- 
trum with uje 



29 



to 



200meV was shown by Park et al 
provide a phonon-induced electron self-energy in excellent 
agreement with that from full first principle calculations. 
Note that we do not take into account corrections to the 
phonon energy and lifetime provided by the phonon self- 
energy as modified by interaction with Dirac electrons. 



These corrections have been calculated theoretically 5|30 
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and their effects have been seen in experiment 3l][32 . For 
the most part, the effects are small for our purposes. Ra- 
man experiments 31 have found a 5meV modulation of 
the energy of the optical phonon line in magnetic fields 
up to 30T. In graphite samples with enhanced lifetimes 
of Dirac fermions, experiments 32 find modulations of 
about 8meV for fields around 5T. Also, the measurements 
see splitting of the phonon structure into three peaks due 
to magnetoresonance of the phonon line energy crossing 
the Landau level. This feature may be relevant for our 
discussion when we examine experimental data, and we 
will return to it then. 

Equations ([I]) and ^ are solved iteratively, beginning 
with the bare density of states, obtained by setting £ — > 
in eq. ([!]). This zeroth-order approximation N^(lo) has 
peaks at w = M n — for all n (up to the cutoff deter- 
mined by W). The dashed red curve in fig. [fja) shows 
this zeroth-order approximation for zero chemical poten- 
tial (the case of charge neutrality) and a f5T field, the 
largest field strength used in tunneling experiments to 
date [7{[8{ [25|j27] . To determine the locations of the peaks 
in the renormalized DOS, we evaluate the imaginary part 
of the first-order approximation immediately found 

to be 
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This carries an image of the DOS, but displaced by loe 
for lo > loe and by —lue for lo < —loe- Therefore, 
the peaks at u = M n — [i in N^(lo) induce peaks at 
(1q ± loe in S^, as seen in the lower frame 



lo 



of fig. [lja). E^(cj) is peaked at these same energies, as 
shown in the middle frame of fig. [fja). The renormal- 
ized DOS iVW(cj) is obtained by substituting EW(cj) into 
eq. ([!]). Since the right-hand side of eq. ([I]) contains not 
just peaks, but a slowly increasing portion contributed by 
the sums of tails of Lorentzians, peaks in S(w) introduce 
new peaks in N(lo) at the same energies. Therefore, after 
renormalization, N^(u) has peaks at lo = M n — fiQ±Lo E , 
which we shall refer to as phonon peaks. This is evident 
in the top frame of fig.[TJa), where we see phonon peaks at 
loe = 200mcV and oje+Mi — 340meV in the renormalized 
DOS (solid black curve). The magnitude of these peaks is 
proportional to the electron-phonon coupling strength A, 
to which we have given the small value of A = 50meV in 
fig. [T] to make for easy comparison between the bare and 
renormalized curves. 

Along with these new peaks, there are those correspond- 
ing to the original Landau levels, now located at solutions 
to lo — Si(w) + /i = M n . The n = LL does not shift, 
remaining pinned at the DP [note that the n = LL 
of the bare DOS is hidden under the black curve in fig. 
1(a)]. For small lo, —Ei(w) ~ lo\ cS , where A cff is the ef- 
fective mass renormalization parameter. As seen in the 
middle frame of fig. [lja), this linear behavior remains a 



good approximation until lo approaches the peak at loe- 
For the parameters used in fig. [lja), A off = 0.04 and the 
n — 1 LL is appropriately shifted down by a factor of 
1 + A off , which can be interpreted as a renormalization of 
the Fermi velocity in M n . But the n — 2 level, which falls 
near loe, is shifted by a larger factor. Between the first 
and second phonon peaks, S(w) returns to linear behav- 
ior. In addition to their locations, the widths of the LLs 
are affected by the self-energy: the imaginary part of £ 
becomes nonzero above loe due to the opening of a new 
scattering channel, broadening the levels above loe beyond 
their bare-band widths. 

In fig. [ljb) , we show a sequence of curves for the renor- 
malized DOS at five increasing values of B between 7.6T 
and 9.12T. The other parameters are identical to those in 
fig. [lja). For clarity each curve is displaced upwards by 
0.35iV o eV. This progression shows what happens when a 
given LL (indicated by the black arrows) is made to pass 
through the phonon energy. In the lower curve (black), 
the LL falls below the small phonon peak that begins 
sharply aX loe- As we increase B, the LL passes through 
the phonon peak. In the process, it is split in two, as in- 
dicated by the two black arrows, one to each side of loe- 
Because the phonon peak is small for our chosen value of 
A, it is overwhelmed by the right peak in the split LL. At 
sufficiently high B, the left peak vanishes and the right 
peak moves sufficiently far from lo e for the phonon peak 
to again become visible. 

We illustrate the effects of renormalization on the en- 
ergy levels in the schematic fig. [ljc) . Comparing the bare 
LLs (shown on the left Dirac cone) to the renormalized en- 
ergy levels (on the right Dirac cone) , we see the LLs (disks 
outlined in dark blue) brought to lower energies, the new 
phonon peaks (dashed red circles) interspersed between 
them, and LLs split in two (here n = ±2 — > ±2 ± ) when 
they are sufficiently close to a phonon peak. 

Our results' experimental relevance is seen by compar- 
ing them with the results of Miller et al. [8], who obtained 
DOS curves from scanning tunneling spectroscopy on epi- 
taxial graphene grown on the carbon-face of SiC. In fig. [2] 
we show a sequence of theoretical DOS curves calculated 
for values of magnetic field in the range to 6T, as labelled 
in the figure. (Curves are again offset for clarity.) Here, 
through consideration of the Miller data, we have used 
parameters A = 0.07, lo e — 230meV, T = 3meV, and 
vf = 1-13 x 10 6 m/s. Note that we have also used a finite, 
though small, chemical potential in order to roughly mimic 
its increase with B seen by Miller et al.; from bottom to 
top, n = 3.4meV, 2.8meV, -3.8meV, -6.6meV, -4.8meV, 
-8.6meV, and -13.3mev. This shift in fi with magnetic 
field has been attributed to a redistribution of charge in 
the graphene multilayers, screening, and issues associated 
with the probe tip [8] , and it results in a visible asymmetry 
between the peak structures at loe and — lo E - Clearly, as 
expected, the effects of phonon coupling grow along with 
the increasing amplitude of the LLs as B increases, while 
for B — (continuous black curve) they cannot be seen on 
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Fig. 2: (Color online) Density of states as a function of energy 
for a sequence of magnetic field strengths. Peaks are labeled 
with their Landau level indices. The dotted lines mark plus and 
minus the phonon energy. N is in units of iVoeV and curves 
are offset by 0.5iVoeV for clarity. 



the scale of the figure. As in the preceding discussion, be- 
low the phonon energy (indicated by the dotted lines), the 
LLs behave as the bare DOS, but renormalized by 1 + A cff . 
Above we, we see the first phonon peak more prominently 
as B increases, followed by more ordinary LLs at higher 
energy, but with heavily damped amplitudes; these effects 
are larger than in the results of Fig. [T] due to the larger 
value of A (or that of X eS ). All these features are in qual- 
itative agreement with the data of Miller et al. , who show 
a similar figure for the same magnetic field values. 

In fig.|3j we show the energies at which the various peaks 
fall (points connected by solid lines to guide the eye) as a 
function of \JB. The parameters are as in the preceding 
figure, but for simplicity, here we use /x = for all values 
of B. For small values of magnetic field, the energy of the 
nth level agrees well with M n /(1 + A cff ) (dashed lines). 
As expected, however, when they approach the phonon 
energy u>e (dotted line), the LLs deviate from this behav- 
ior and are split in two about uje- In the inset, we show 
that when plotted as a function of %/ nB, all the curves col- 
lapse into a single straight line of slope v p y // 2he/c, except 
for the notable deviation near the phonon energy. 

Figure [4] shows a more detailed comparison with the 
data (black dots) of Miller et al. [8j in the case B — 5T. 
The solid blue curve is our result for the renormalized DOS 
(in arbitrary units) as a function of energy w, still with 
the same parameters as in fig. [2| We see qualitative and 
even semiquantitative agreement (left frame). The phonon 
structure around 230meV is clearly seen in the data, as is 
a clear drop in amplitude of the LLs beyond this energy. 
In the right frame, we show an enlarged view of the re- 
gion around the phonon energy. Along with results for 
coupling to a single Einstein mode (dashed red curve) , we 
show results for the more realistic model of a truncated 
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Fig. 3: (Color online) Locations of peaks in the DOS as a 
function of yB. The circular data points connected by solid 
curves show the positions of the renormalized Landau levels. 
The thin dashed lines show M n /(1 + X cS ). The dotted black 
line indicates the phonon energy. Curves for negative energies 
are the mirror image of those shown here, fnset: the same 
points plotted as a function of V nB. 



Lorentzian phonon distribution, centered at 230meV with 
a half-width of lOmeV and a cutoff at \u> — oje\ — 20meV. 
In the former case, the n — 7 Landau level is sharply split 
around the phonon peak; in the latter case, the phonon 
peak is broadened and merges with the split Landau level, 
creating the flattened structure seen in the data. This 
structure lies at an energy roughly 30meV greater than 
that of the i?2g mode. A phonon structure at a similar 
energy, 240meV, was also identified by Brar et al. [12] in 
their STS study of graphene flakes on Si02 with no mag- 
netic field, and they suggest that it could be due to a 
multi-phonon process. Alternatively, it might be the E 2g 
mode shifted to higher frequency by some other mecha- 
nism. For example, tip effects such as deformation of the 

and other anomalous fea- 



sample could cause a shift 33 



tures in Miller's data have been ascribed to them [34] . A 
shift could also be caused by the magnetoresonance ef- 



fect seen in Raman measurements 31 32 , which is espe 



cially large at 5T. Regardless of the origin of this feature, 
our model provides a reasonable fit to the data and gives 
a rather small value of effective mass renormalization of 
A cff ~ 0.07. 

In conclusion, we have calculated the effect of electron- 
phonon interactions on the density of states N(lo) in 
graphene in a magnetic field B. The imaginary part of 
the electron-phonon self-energy gives the electron scatter- 
ing rate associated with an electron scattering between 
states of energy uj and w-u^, where uje is the phonon 
energy. This process is possible only for lo > lue, and it 



p-4 



Phonon structures in the electronic density of states of graphcne in magnetic field 




I £ I i I i I I I i I i I i I 

100 200 300 200 220 240 260 



w(meV) Lu(meV) 

Fig. 4: (Color online) Density of states as a function of energy 
at B = 5T, comparing numerical results to experimental data 
(solid black circles) taken from Miller et al. [8]. The plot on 
the right zooms in on the region around the phonon energy, 
showing results for both the Einstein model (dashed red curve) 
and a truncated Lorentzian (solid blue curve) . 



carries an electronic density of states factor N(ui— u>e)> en- 
coding information on the underlying LL structure. From 
this encoding, the electron-phonon coupling generates a 
new set of peaks, each one an image of a bare Lan- 
dau level displaced by the phonon energy. As the mag- 
netic field is increased, both the LL peaks and the new 
peaks are increased in amplitude. At energies small com- 
pared to W£, the renormalized Landau levels closely fol- 
low E n = sgn(n)y/2\n\eBv 2 F h/c/(l + A cff ), where A off is 
the electron-phonon mass renormalization parameter. As 
uje is approached, however, the LL is shifted beyond the 
1/(1 + A eff ) renormalization factor; furthermore, it is split 
in two about the new peak at u>e- Since for uj > lue a new 
scattering channel opens, the resulting amplitudes of the 
LL lines above loe are also reduced over their bare band 
value and their widths increased. All these features are 
seen in the STS data. Fitting to this data at B — 5T 
provides a new estimate of the electron-phonon mass en- 
hancement A cff ~ 0.07. 
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